Transition from isolated to overlapping resonances 
in the open system of interacting fermions 
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We study the statistical properties of resonance widths and spacings in an open system of inter- 
acting fermions. At the transition between isolated and overlapping resonances, a radical change in 
the width distribution occurs with segregation of broad ( "super-radiant" ) and narrow ( "trapped" ) 
states. Our main interest is to reveal how this transition is influenced by the onset of chaos in the 
internal dynamics regulated by the strength of random two-body interaction. In the transitional 
region, the width distribution and its variance, as well as the distribution of spacings between reso- 
nances are strongly affected by internal chaos. The results may be applied to the analysis of neutron 
cross sections, as well as in the physics of mesoscopic devices with strongly interacting electrons. 
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Physics of marginally stable mesoscopic systems is an 
important part of current scientific interest, both in nu- 
clear physics and in condensed matter physics. The cen- 
ter of attention in nuclear physics moved to nuclei far 
from the valley of stability. Exotic nuclei have low bind- 
ing energy so that the continuum states are easily excited 
by a weak external perturbation and appear as the closest 
virtual states in quantum- mechanical calculations. The 
correct description of interplay between reactions and in- 
trinsic structure is also important for practical applica- 
tions connected to neutron cross sections at low and in- 
termediate energy. Mesoscopic solid-state devices, such 
as quantum dots and quantum wires, are open being ex- 
plicitly coupled to the outside world. The problem of 
mesoscopic conductance fluctuations |1] as a part of more 
general physics of quantum transport in many aspects is 
similar to that of fluctuations in nuclear reactions Q. 
In all such transition from a consideration of 

a closed model system to a realistic open system in its 
interaction with continuum becomes necessary. 

A working instrument in the unified description of in- 
trinsic states and reactions, or decay channels, is an ef- 
fective non-Hcrmitian Hamiltonian TL derived by elimi- 
nation of channel variables, for example with the help 
of the projection method [1, 0> S B 0- This Hamil- 
tonian describes the intrinsic dynamics of an open sys- 
tem; its eigenvalues are complex energies, Sj = Ej — i Tj , 
where the width T determines the lifetime of a resonance, 
r ~ h/T. The form of the imaginary part of Ti, factor- 
ized in the amplitudes of entrance and exit channels, is 
dictated by the unitarity of the scattering matrix. The 
factorization brings in remarkable consequences observed 
in numerical simulations [8] and explained theoretically 
in Refs. where the analogy to optical superradiance 
[irj| of coherently radiating atomic emitters was pointed 



out. 

When typical widths T are small compared to the level 
spacing, D, on the real energy axis, T/D <C 1, the cross 
sections display isolated resonances. In the critical region 
T/D « 1 with crossover to overlapping resonances, the 
width distribution displays sharp segregation of broad 
short-lived (superradiant) states and very narrow long- 
lived (trapped) states @, [tl [ll|. Correspondingly, the 
distribution of poles of the scattering matrix undergoes 
a sharp transition from one to two "clouds" of poles 
in the complex plane of resonance energies 12]. This 
crossover phenomenon is universal. However, it is not 
yet known how the specific features of this universal ef- 
fect, and therefore the characteristic time scales, depend 
on the internal interactions. 

The region of overlapping resonances was and remains 
of special interest starting from the seminal papers by 
Ericson [l3| where the general approach to fluctuations 
of cross sections has been developed. In nuclear physics, 
this region until now was not a subject of detailed study, 
partly because of lack of a statistically reliable experi- 
mental body of data, however, it is accessible in nuclear 
reactions for scattering of slow neutrons on heavy nu- 
clei. In solid states physics, strongly overlapping reso- 
nances occur in the regime of the so-called perfect cou- 
pling, where famous universal conductance fluctuations 
arise [l|. The similarity and differences between con- 
ductance fluctuations and Ericson fluctuations in nuclear 
reactions have been discussed in Ref. 0. Recently, ex- 
periments in the region of overlapping resonances have 
been performed in microwave cavities, see Ref. [Til. [l5j . 

Below we study the statistical properties of resonance 
widths and spacings in the transitional regime explic- 
itly taking into account the regular or chaotic dynamics 
inside the system [l|| Il7j . The local statistical prop- 
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erties of fully chaotic quantum states can be identified 
with those of the Gaussian Orthogonal Ensemble (GOE) 
fl8j |. Here, for the first time, the two-body random in- 
trinsic interaction is considered that allows one to vary 
the degree of internal chaos in an open system. With 
the non-Hermitian ap pro ach, random two-body dynam- 
ics was used in Ref. [13 ] in the limit of a very strong 
interaction, when the influence of the regular mean field 
was neglected. 

We consider a set of a large number, N, of intrin- 
sic many-body states \i) with the same exact quantum 
numbers. The states are unstable being coupled to M 
open decay channels. The dynamics of the whole system 
is governed by an effective non-Hcrmitian Hamiltonian 
[Hill given by a sum of two N x N matrices, 

M 

H = H- l -W; I^=$>^ c ; (1) 

the intrinsic states are labeled as i, j, ... and decay chan- 
nels as a,b,c... Here H describes Hermitian internal dy- 
namics, while W contains the amplitudes A\ coupling 
intrinsic states \t) to the open channels c. Under time- 
reversal invariance, both H and W are real symmetric 
matrices. 

The Hermitian part, H = Ho + V, of the full Hamilto- 
nian |T]) is modeled by the so-called two-body random in- 
teraction (TBRI) of n fermions distributed over m single- 
particle states; the total number of many-body states is 
N = m\/[n\(m — n)!]; in our simulations n = 6, m = 
12, N — 924. In this model Ho describes the mean 
field part, where single-particle energies, e s , are assumed 
to have a Poissonian distribution of spacings, with the 
mean level density 1 /dp . The two-body interaction V 
between the particles [17| is fixed by the variance of the 
two-body random matrix elements, (V^ S2; s 3 ,s 4 } = v o- 
While at vo = we have a Poissonian spacing distribu- 
tion P(s) of many-body states, for do = (infinitely 
strong interaction, vo/do — > oo), P(s) is close to the 
Wigner-Dyson (WD) distribution typical for a chaotic 
system. The critical interaction for the onset of strong 
chaos is given [llj by v cr /do ~ 2(m — n)/N s , where 
N s = n(m — n)[l + (n — l)(m — rt — l)/4] is the number of 
directly coupled many-body states in any row of the ma- 
trix Hij. Thus, we have in our model v cr /do ~ 1/20, and 
often we have used the value vo/do — 1/30 slightly less 
than v cr . In parallel we also consider the limiting case 
of H as a member of the GOE [2(| that corresponds 
to a many-body interaction, when the matrix elements 
are Gaussian random variables, (Hfj) = l/N for i ^ j 
and (H?j) = 2/N for i = j. In order to compare with 
the GOE case, when changing the strength vq of interac- 
tion, we renormalize the model parameters to have the 
same mean level spacing between many-body states at 
the center of the resulting energy spectrum. 

The real amplitudes A\ are assumed to be random in- 
dependent Gaussian variables with zero mean and cor- 
relator l \A\A C j ) = 6ijS cc -f c /N. This is compatible with 



the GOE or TBRI model assumed for intrinsic dynam- 
ics when generic many-body states coupled to contin- 
uum have complicated structure, while the decay probes 
specific simple components of these states related to few 
open channels. Even in the case of weak intrinsic inter- 
action, we need to have in mind that the states \i) have 
certain values of exact constants of motion, such as angu- 
lar momentum and isospin in the nuclear case. At suffi- 
ciently lar ge d imension N, these states acquire geometric 
chaoticity [211 ] due to the almost random coupling of indi- 
vidual spins. Therefore the ensemble of decay amplitudes 
is reasonable. The parameters j c /N with dimension of 
energy are unperturbed partial widths which characterize 
the coupling to the channel c. The normalization of this 
correlator is convenient if the energy interval ND covered 
by decaying states is finite. We neglect a possible explicit 
energy dependence of amplitudes (the transitional region 
is far from thresholds). 

The effective Hamiltonian (JTJ) determines the reaction 
cross sections, a ba (E) oc |>S ba (-E)| 2 , or the scattering ma- 
trix S = 6 ba — iT ba . The reaction amplitudes, neglecting 
the smooth potential phases irrelevant for our purposes, 
are expressed in terms of the amplitudes A\: 

N ( 1 \ 
T b *(E)=Y^AU— -J A*. (2) 

The complex eigenvalues £j of TL coincide with the poles 
of the S'-matrix and, for small 7 C , determine energies 
and widths of separated resonances. In the simulations 
we consider an energy interval at the center of the spec- 
trum of H. with the constant many-body level density 
p(0) = D^ 1 . As 7 C grows, the resonances start to over- 
lap, the effective parameter being k c = irj c /2ND. The 
transmission coefficient in the channel c, 

A K C 

T< = 1-|<^)| 2 = ITW (3) 

is maximal (equal to 1) at the critical point, k c — 1, that 
marks the transition to superradiance and trapping. We 
study the statistical properties of resonance widths and 
spacings as a function of the interaction between particles 
(ratio vo/do) and the continuum coupling parameter k. 
For simplicity, we assume M equiprobable channels, k c = 
k; the maximum value of M we considered was M = 10. 
For each value of k we have used N r = 100 realizations 
of the Hamiltonian matrices, with further averaging over 
energy. 

At a critical value, k » L a segregation of the reso- 
nance widths occurs 

H S El, see inset in Fig. [Tj The 
widths of M resonances are increasing at the expense of 
the remaining N — M resonances. For weak coupling, 
k < 1, the widths are given by diagonal matrix ele- 
ments, Ti — (i\W\i) — J2cLi(Ai) 2 , and the mean width 
is (r) = jM/N. In the limit of strong coupling, k» 1, 
the widths of M broad resonances converge to the non- 
zero eigenvalues of the matrix W that has a rank M due 
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to its factorized structure which is dictated by unitarity 
of the scattering matrix. As for the remaining "trapped" 
(N — M) states, their widths decrease oc Ih. 




FIG. 1: (Color online) Average width versus the coupling 
strength n for M = 10. Solid curves show the expression 
@, pluses refer to vo = 0, crosses to vq = do/30, squares to 
vo/do = oo. In all cases for k > 1 the average was made over 
N — M narrow resonances only. More details are shown for 
vo = do/30 in the inset, where open circles correspond to the 
average over all N resonances, and full circles stand for the 
average over all resonances for k < 1, and over N — M narrow 
resonances for k > 1. 

In Fig. [T] we show how the average width (normalized 
to the mean level spacing D) depends on k for different 
values of vo/do. The inset demonstrates a sharp change in 
the width distribution at the transition point K = 1. The 
results have been compared with the Moldauer-Simonius 
(MS) relation [22], |23j (that follows also from our statis- 
tical assumptions) 
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For our case of M equivalent channels, 

M - — 1 + K 

D ~ tt n 



1 - K 



(4) 



(5) 



Here we put the absolute value of the ratio under the 
logarithm in order to extrapolate this expression beyond 
the transition point, for k > 1. The extrapolated MS- 
formula works well for any interaction stren gth . For one 
channel, this result was obtained in Refs. [jjl Till]; nere we 
see that the MS-expression is also valid for a large number 
of channels, independently of the interaction strength. 
This agrees with the well known fact that average quan- 
tities are in fact functions of the transmission coefficient 
that is directly related to k. This statement is obvious 
for the GOE intrinsic dynamics where the results have 
to be expressed in terms of orthogonal invariants. In our 



case this follows from the assumed statistics of the con- 
tinuum amplitudes when we make this averaging first, 
independently of the degree of intrinsic chaoticity [H| • 
According to [HI], the divergence of (r) at n = 1 is due to 
the (non-integrable) power-law behavior for large T, see 
below; in the numerical simulation, see Fig. [TJ there is 
no divergence because of the finite number of resonances 
although the trend is clearly seen. 
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FIG. 2: (Color online) Width distribution for M = 10 and 
k = 0.01; 1; 10, with the same symbols as in Fig. [1] Smooth 
curves are the Xio distribution. 

The width distribution P(T) is shown in Fig. [2] for dif- 
ferent values of vo/do and k = 0.01, 1, 10. At weak cou- 
pling, Fig. 2a, the conventional ^/-distribution is valid 
for any strength of the interaction. However, when k 
increases, Fig. 2 b and 2c, a clear dependence on the in- 
teraction strength emerges. As first noted in Ref. (26|, as 
k increases, P(T) becomes broader than the Xm distribu- 
tion. For large k, both for the GOE and for vo/do — > oo, 
P(r) is again given by the Xm distribution, contrary to 
the cases of the finite interaction strength; this reflects 
the uniformity of properties of all fully chaotic intrinsic 
states. 

Apart from the one-channel case studied in Ref. a 
general analytical treatment of the width distribution is 
a difficult task even for simplified random matrix mod- 
els. In the case of GOE for the Hermitian part H, the 
width distribution was found in Ref. [25[, neglecting the 
emergence of broad resonances in the regime of strong 
overlap at k > 1. Our results for vo/do = oo correspond 
well to the analytical predictions. In particular, the tails 
of -P(r) decrease ~ T~ 2 that leads to the divergence of 
(r) in agreement with the MS-formula. 

The (normalized) variance of the resonance widths, see 
Fig. 3 for M = 10 and k = 0.8, strongly depends on the 
interaction between particles. Our results confirm the an- 
alytical estimate for the transition to strong chaos that 
can be associated with the WD-distribution of level spac- 
ings. Note, however, that the transition is quite smooth, 
and even for a relatively strong interaction there is a de- 
viation of the variance from the GOE value. We would 
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like to stress that the level spacing distribution, as the 
weakest signature of quantum chaos [16], turns out to 
be rather insensitive to the transition. Indeed, the data 
in the inset for n — roughly correspond to the WD- 
distribution, although once the continuum coupling is 
switched on, the variance of the resonance widths is very 



is the width variance. For the GOE and vo/do — > oo, the 
system returns to the \m distribution, with an increase 
of k. This is not true for a finite value of vo/do, in agree- 
ment with the results of Ref. (28|, where justification of 
the broadening of the width distribution in the Poisson 
case (equivalent to vq = 0) as compared to the GOE was 
given for M = 1. 




FIG. 3: (Color online) Dimensionless variance of the widths 
versus vo/do for M = 10 (connected circles): the GOE value 
is shown by horizontal dashed line; the vertical dot-dashed 
line marks the value vo = v cr corresponding to the onset of 
chaos. In the inset the level spacing distribution is shown for 
ft = and vo/do = 0.2 (crosses), see the arrow in the main 
part, and for the GOE (circles). The smooth curve is the 
WD-distribution. 
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FIG. 4: (Color online) Dimensionless variance of the width vs 
k for M = 10. The dependence due to the xl/ distribution is 
shown by the horizontal line; the symbols are the same as in 
Fig. [1] Non-smooth dependence for n > 1 is due to a shift of 
the critical point for M » 1 from n = 1 to K = K c > 1. sec 
Ref. Q. 

It is instructive to show how the variance of the widths 
depends on the continuum coupling, Fig. [4] The pre- 
diction of the Xm distribution that the variance is 2/M 
independently of k is correct for very weak and very 
strong coupling. As expected, the variance takes its max- 
imal value at the transition point n — 1 (it diverges for 
N — > oo). As for the dependence on vq, the less chaotic 
is the intrinsic dynamics and therefore greater a possible 
difference in the structure of decaying states, the larger 




0.01 0.1 1 10 




0.01 0.1 1 10 



FIG. 5: (Color online) The probability P(0) of finding the 
level spacing s < 0.04 (between the real parts of poles of 
the scattering matrix, in units of mean level spacing) vs k. 
Upper panel: GOE for M = 1 (open circles) and M = 10 
(full circles); lower panel: M = 10 with the same symbols as 
in Fig. [T] The value due to the WD-distribution is given by 
the horizontal line. 

Finally, we analyze how the level repulsion (in our case 
the repulsion of the centroids of resonances along the real 
energy axis) depends on the continuum coupling and on 
the intrinsic interaction. The characteristic quantity is 
the probability of the level proximity, P(s — ► 0). It is 
finite for the Poisson distribution and vanishes linearly 
for the GOE. A weakening of the level repulsion with 
an increase of the continuum coupling (or the number of 
channels M) was noted as a generic property in Ref. . 
The physical explanation (energy uncertainty of quasis- 
tationary states) was given in |9j. Our detailed data, 
Fig. O confirm these results and show the dependence on 
vo/do. For non- interacting particles, vq = 0, the results 
are consistent with those of Ref. [29( , where it was argued 
that, for the Poissonian intrinsic Hamiltonian H , the dis- 
tribution P(s) becomes similar to that of a chaotic sys- 
tem due to the perturbation induced by the coupling to 
the continuum. Our data clearly show that with increase 
of intrinsic chaos the level repulsion grows, however, for 
finite i>o always remains weaker than for the GOE case. 

In conclusion, we have studied statistics of complex 
energies (resonance widths and spacings) for a generic 
fcrmion system coupled to open decay channels. For the 
first time we carefully followed various signatures of the 
crossover from isolated to overlapping resonances in de- 
pendence on the inter-particle interaction. We show that 
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the Moldauer-Simonius expression for average widths is 
valid also in the region k > 1, provided the M widest res- 
onances are excluded. Also, we found that the average 
widths are insensitive to the interaction strength. On the 
other hand, the width variance depends on the degree of 
intrinsic chaos at the transition point /t » 1 in a non- 
trivial way. Even if for a closed system (k = 0) the level 
spacing distribution is of the Wigner-Dyson form, the 
width distribution for a perfect coupling to the contin- 
uum (k « 1) is very different from the GOE predictions. 

Our study is based on the two-body random nature of 
the inter-particle interaction. However, since we explored 
the central part of the energy band where the eigenstates 
are delocalized, one can expect similar statistical prop- 
erties of resonances if one substitutes the intrinsic inter- 
action V by the GOE matrix, with a proper rescaling of 
the interacting strength vq. This may allow for an accu- 
rate analytical analysis. Another extension should con- 
sider realistic shell-model interactions which bring, at a 
high level density, the results close to those for banded 
random matrices. The results on statistical properties 
of the reaction cross sections, including deviations from 
the theory of Eriscon fluctuations, are reported elsewhere 
[24j . The obtained results can be applied to neutron res- 
onances in nuclei and to open mesoscopic systems in the 



crossover region, where we found strong deviations from 
the conventional level statistics and \ 2 width distribu- 
tion. Although it might be difficult to extract from the 
experiment the exact statistics of poles in the complex 
plane (which would be extremely desirable), at least the 
trends of distributions along the road from isolated to 
overlapping resonances certainly can be studied in well 
measured neutron resonances on heavy nuclei at energies 
further away from threshold. The experiments with mi- 
crowave billiards (analog of one-body shape-dependent 
chaos) at variable continuum coupling can shed light on 
some features of the process of width collectivization and 
segregation. Similar effects are expected in mesoscopic 
conductance fluctuations which will be discussed in our 
forthcoming publication. 
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